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Abstract 

We review the current status of the problem of chiral symmetry restoration, discussing 
an open question of whether at the SU{Nf)A restoration point T = Tc the U{1)a chiral 
symmetry is or is not (approximately) restored. New lattice and instanton-based studies 
are considered, and some new calculations are suggested to clarify the issue. We also 
speculate on possible experimental manifestations of two possible scenarios. 



1 



l.The problem 

For simplicity, we ignore all effects due to the non-zero quark masses, and consider 
QCD in the chiral limit, with Nf massless quarks. In this case the QCD Lagrangian is 
just a sum of two separate terms, including right- and left-handed quarks, which implies 
two chiral symmetries: SU^Nf )^ and U{1)a- 

Their fate is well known to be different. The former one is spontaneously broken 
in the QCD vacuum but it is restored at high temperatures, above some critical point, 
denoted as T = Tc. Present lattice simulations with dynamical quarks suggest ~ 
ISOMel/. 

The U{1)a chiral symmetry is not related to Goldstone bosons, as Weinberg has 
first pointed out. It is now well known this is is because this symmetry simply does not 
exist at quantum level, being violated by the 'chiral anomaly'. Its physical mechanism 
is also known it is driven by the tunneling between the topologically different gauge 
vacua, described semiclassically by the instantons 0. 

It is also known, that at high temperatures the instanton-induced amplitudes are 
suppressed due to the Debye-type screening and therefore (at some accuracy level) 

we expect this symmetry to be 'practically restored' at high T. Let us denote the point 
where it happens with some reasonable accuracy as 

The main question to be discussed below is the interrelation of the two temper- 
atures, Tc and Tu{i). Let us refer as 'scenario 1' to the case Tc -C in which the 
complete U{Nf)A chiral symmetry is restored only well inside the quark-gluon plasma 
domain. Another possible case[| ~ ^c/(i) is referred below as 'scenario 2': it implies 
significant changes in many hadronic channels around this phase transition point. As we 
will discuss below, these two scenarios lead to quite different predictions. 

This important question is well known and it was already discussed in literature, 
but still we do not have a definite answer. In these comment, we are going to look again 
at its somewhat controversial history, and discuss what was recently done (and also can 
be done) in theory, numerical and 'real' experiments in order to answer it. 

^The case Tc >> does not seem to be possible. 

2 



Pisarski and Wilczek ||^ have considered this question in connection with the 
order of the chiral phase transition. They have pointed out that in the special case 
Nf = 2 the 'scenario 1' is hkely to lead to the second or(ier transition. The reason is an 
effective Lagrangian describing the softest modes is essentially the Gell-Mann-Levy sigma 
model, same as for the 0(4) spin systems. 

However, in that paper Pisarski and Wilczek have actually argued in favor of the 
'scenario 2'. Their argument was as follows: 'if instantons themselves are the primary 
chiral-symmetry-breaking mechanism, then it is very difficult to imagine the unsuppressed 
t/(l)A-breaking amphtude at Tc'. They have even mentioned that this amplitude should 
be at Tc at least an order of magnitude smaller than at T=0, although no details of this 
estimate were given. 

As we will discuss below, both lines of reasoning have accumulated new evidences, 
so the question under consideration appears to be even more controversial, than it was a 
decade ago. We hope new studies, especially the lattice-based ones, can finally clarify it. 

2.1s the Nf = 2 phase transition similar to that for 0(4) spin system? 

During the last 2 years the large-scale numerical lattice simulations with dynam- 
ical fermions have addressed the issue of QCD phase transition, especially its order. It 
was indeed found that the case of 'many fiavors' Nf > 2 (especially well studied case 
is 4) leads to sharp first order transition, with two phase coexistence etc., while in the 
particular case Nf = 2 such phenomena were not observed (see discussion and original 
references in |^ or proceedings of yearly lattice conferences). 

Rajagopal and Wilczek [0] have argued that in this case one probably has a second 
order transition, analogous to that of the 0(4) spin system. They have therefore proposed 
a suitable effective Lagrangian, essentially a 3-dimensional a model, and have applied it 
to dynamical simulations of the phase transition. Now, do the available lattice data really 
support these ideas? 

The most straightforward way to test them is to compare the critical behaviour 
in both cases, testing whether the Nf = 2 QCD and the 0(4) spin system do or do not 
belong to the same universality class. (Similar comparison was done in the past for many 
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different phase transition, including 'deconfinement' one in SU{2)c pure gauge theory.) 

The first critical index to compare is the one for the order parameter, for which 
the analogy ||^ suggests 

<^^>^|(T-T,)/T,|-" (1) 

Unfortunately, it is not that simple to compare it to lattice data, because chiral order 
parameter is only obtained as a small- mass limit, and the masses used on the lattice 
cannot be made sufficiently small. Recent analysis has concluded, that the data are 
consistent with 0(4) critical exponents, although say 0(2) ones are not also excluded. 

The second obvious question is the behaviour of global thermodynamical quanti- 
ties, such as the specific heat. The 0(4) spin system is believed to have a very amusing 
behaviour, with positive power^ 

C(T)~|(T-T,)/T,|-^9±-06 (2) 

It means that the singular contribution of the soft modes vanishes at the critical point, and 
in order to single it out the 3-ed derivative of the free energy should then be calculated. 

However, lattice data for the Nf = 2 QCD actually do show a huge peak in the 
specific heat around Tc. It certainly implies, that many new degrees of freedom become 
available (or are significantly changed) in this region. What these degrees of freedom are, 
both in hadronic language and in the quark-gluon one, remains the major open problem 
in the field. 

Of course, there is no logical contradiction here: apart of large (but smooth) peak 
one may eventually find a small 'kink', which is truly singular. And still, if one can get 
any hint from the available lattice data at all, I think they certainly indicate that point 
to the o the Nf QCD and the 0(4) spin system are very different physical objects, even 
as far as such global parameters as thermodynamical properties are considered. 

3. Are both chiral symmetries restored simultaneously? 

^As far as I know, it remains unknown whether the coefficient is positive or negative: thus one can 
have a dip or a peak. 
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We have already mentioned at the end of section 1 the arguments in favor of 
'scenario 2' mentioned in ||^. Are there any new arguments for or against it? 

The early conjuncture (see e.g.[Q and references therein) relating instantons with 
spontaneous breaking of the SU{Nf)A chiral symmetry has gained strong support during 
the last decade. From a qualitative idea one has come all the way to quantitative under- 
standing of many details. It is now possible to to calculate not only the quark condensate, 
but also such quantities as the pion or nucleon masses, using instanton-based model. 

The latest development has resulted in quite detailed studies of the QCD corre- 
lation functions. The 'instanton liquid' model was shown to reproduce their behaviour 
for many mesonic and baryonic channels in agreement with phenomenology JlOl and 



lattice calculations |TT|. Direct lattice studies [|T^, have also shown direct connection 
between the quark condensate (and correlation functions) and instantons. 

Thus, there is little doubt that instantons drive chiral symmetry breaking in the 
QCD vacuum. Unfortunately, all this development corresponds to T=0, while very little 
was done for the finite T case. 

Although the Debye screening of large-size instantons 0] should strongly suppress 
them at high T, what exactly happens below Tc remains unknown. The so called 'Pisarski- 
Yaffe suppression factor' 

dmnstiT) ~ dn,nst{T = 0)exp[-7rVT2(^ + ^)] (3) 

can be applied only in the quark-gluon plasma phase^. Taking into account available 
lattice data on thermodynamics and Debye screening length, one can now argue that 
this relation is applicable above T ^ 2Tc. We also know now that the typical instanton 



radius of the 'instanton liquid' is p ~ 1/3 fm |T^. Combining the two one finds instanton 
suppression by at least 2 orders of magnitude at T ^ 2Tc. 

Now, is it possible that strong (e.g. by one order of magnitude, as in the original 
Pisarski-Wilczek estimates) instanton suppression is already present at T = Tc"? There 
^ Whatever simple this argument is, it was not made in the original papers. For small size instantons 
't Hooft determinant is calculated correctly even at T=0, but the temperature-dependent corrections to 
it are calculated under specific assumption that the heat bath is nothing but an ideal gas of quarks and 
gluons. Another way to see this is to note that the Debye screening is absent below Tc- 
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are several reasons to think it is not the case. Unfortunately, none of them is direct and 
sufficiently strong at the moment. Let me mention two of them. 

The first are measurements of T- dependence of the topological susceptibility x in 
pure gauge theories. Let me remind, that in absence of fermions, it is more or less direct 



measure of the instanton density. What is found [|T6| demonstrate that x{T) remains 
constant (or even somewhat increases) through the deconfinement phase transition, while 
the high-T data are so far inconclusive. Note however, that deconfinement in pure gauge 
theory corresponds to relatively high temperatures, about 240 MeV, compared to chiral 
restoration we discuss in the rest of the paper. 

The second argument is the following one. The pressure of quark-gluon plasma 
can be written as, p = pQ'>^<^^''s,giuons^rp^ _ B{T), where the first term is perturbative black- 
body type contribution, and the second 'bag term' describes the difference between the 
non-perturbative energy at zero and non-zero T. Strong suppression of instantons at 
implies that the second term is mostly the vacuum energy 

B = 6/1287r2 < {gC^J'^ >^ AOOMeV/fm^ (4) 

(where b is the famous coefficient of the Gell-Mann-Low beta function). Comparing it 
with the perturbative term, one then finds that it is too small to get positive pressure of 
the plasma phase. More details on lattice data and 'melting' of quark condensate can be 
found in ref. where it is concluded that instanton suppression factor at is about 
1/2 or so. Although these arguments are not of course strict, they do indicate that strong 
instanton suppression at Tc is very unlikely. 

Moreover, the qualitative picture of instanton-driven chiral symmetry restoration 
has significantly changed since the days of original conjectures. In fact, suppression of 
instantons is not the only way to 'kill' the quark condensate. Not only the number of 
instantons is important, but also their relative positions and orientations. If the tunneling 
events, instantons and anti-instantons, are rearranged into some finite clusters with zero 
topological charge, such as well- formed 'instanton-anti- instanton molecules', all properties 
of the vacuum are completely changed. In particular, such ensemble leads to zero quark 
condensate, or restored chiral symmetry. 
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Numerical simulations of interacting 'instanton liquid' at finite T shows that 
it is exactly this type of correlations which are rapidly build up in the critical region. 
'Pairing' of instantons allow them in this case to interact all the way around the Matsubara 
torus. This is especially enhanced in case of many light quark flavors, since light quark 
do not propagate well into space direction at high temperatures. 

To work out a detailed theory of such correlations will take some time, so a 
schematic model for this phenomenon was recently proposed and studied in ref. . The 
instanton ensemble is assumed to be a mixture of two components; (i) a random 'liquid' 
and (ii) // molecules Q. It is easy to generate such ensemble, for various values of the 
'molecule fraction' fm = Miii^ksMies ^ ^j-^g main parameter of the model. Behaviour of the 
main order parameter, < gg > on fm is very similar as its T-dependence, measured on 
the lattice. 

It was further found, that different correlation functions depend on it quite dif- 
ferently. Some correlators are rather insensitive to it. For example, the vector one does 
not change by more than about 10 % for all compositions, from = — 1. Even the pion 
correlator shows remarkable stability for /m = — 0.8, with subsequent rapid drop if its 
coupling constant toward fm= 1, at which point it coinsides with its scalar 'relatives' (see 
below). At the same time, some correlators show dramatic sensitivity to fm, especially 
the scalar ones. 

There is no place here for further details, let me now point out the main physics 
of the phenomenon. Its one side was repeatedly emphasized before ( see e.g. Ijl^): as the 
quark condensate 'melts down', so does a 'constituent quark mass' (whatever it means) 
and hadrons in average should become lighter. However, there is another side of the story: 
the instanton-induced interaction between quarks 'melts' as well. In the channels where 
this interaction is attractive, such as p, vr, the two effects work in the opposite directions 
and can cancel each other. However, in the channels where it is repulsive, such as axial 
and 1=1 scalar ones, they work together and lead to much stronger effects. 

4. What lattice observables can further clarify the problem? 



^Note the difference with earher treatment of the problem, e.g in ref. |l8[ : now these two components 
are not treated as two phases, for T < they are present simulteneously . 



Let me start with the most obvious ones. For example, why not measure the 
T- dependence of the instanton density, either 'globally', by the topological susceptibility, 
or 'locally, by 'cooling' or by the fermion method? 

We have already mentioned some results on T-dependence of the topological 
susceptibility in pure gauge theory: those are not yet very accurate, and can be much 
improved. The same question in theories with quark was studied by Bernard et al [pO[| , 



but in this case xiT) is simply proportional to the quark condensate, as Word identities 
imply. 

It is in principle possible to see 'molecular' contribution, substituting 4-volume 
by its surface in the definition of the topological susceptibility 

^(surface H'^surface^oo ^ Q ^ (5) 

Unfortunately, other methods are not very sensitive to 'molecules'. 'Cooling' can 
lead to 'annihilation' of a correlated pair. 'Fermion' method can help only if the pair is 
not too close, so that they still lead to pretty small eigenvalues of the Dirac operator. 

Another question, which was often asked: why lattice people do not measure 
directly the 1]' mass, which is obviously the one most intimately related with the U{1)a 
chiral symmetry and its (practical) restoration. 

The answer^ is technical: in such case it is clear that one should calculate not 
only the so called 'one-loop diagram' 

K,{x,y) = -< Tr[TS{x,y)TS{y,x)] > (6) 

but that of the 'two-loop' one 

K^ix, y) =< Tr[VS{x, x)]Tr[VS{y, y)] >) (7) 

as well, and the latter is much more difficult to evaluate because one has to invert the 
Dirac matrix many times. 

However, one can come around this problem, by looking into the scalar cor- 
relators. Let us discuss for simplicity only two light flavors and use the old-fashioned 

^It seems a strong argument, forbidding measurements of the rj' mass: as for the correlator at not-so- 
large distances, I think the needed statistical accuracy may still be achieved. 
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notations, calling the isoscalar 1=0 scalar channel a a one, and isovector 1 = 1 scalar 
channel a 5 one 0. Under SU{2)a transformations, a is mixed with vr, thus restoration of 
this symmetry at T^. require identical correlators for these two channels. 

Another chiral multiplet is S,rinon- strange, where the last channel is the SU(2) 
version of f]': at T=0 those are very heavy and are not considered in chiral Lagrangians, 
or course. 

On the contrary, the U{1)a transformations mix e.g. tt, 6 type states, and thus its 
'practical restoration' should imply that such type of correlators should become similar. 
Finally, if both chiral symmetries are restored, a simpler statement follows: left-handed 
quarks never become right-handed, therefore all vr, rjnon- strange, o", ^ correlators should be- 
come the same. 

Now it is clear what one should do: to measure the 6 correlator, for which one 
does not need the double-loop diagram [], as a function of T. Comparison with the pion 
correlator will tell us where they become close enough, which defines the under 
consideration. 

Unfortunately, strong confusion was going on in lattice measurements of scalar 
correlators. A number of such measurements was reported |2^, with Kogut-Susskind 



fermions. In this formulation different flavors live on different lattice sites, so the operator 
is quite different for a, 6. Although a operator (and a label) was used, only the one-loop 
part K^{x,y) was calculated. It was somehow implicitly assumed, that the two- loop one, 
which leads to separable contribution 

K2ix, y) ^< qq{x) >< qq{y) > (8) 

at large distances \x — y\, is something like a constant, which should be subtracted. 

However, there is no ground for such an assumption, and experience with a 
correlator calculated in the 'instanton liquid model' P] shows that, even after subtraction, 
there remains a significant contribution of some light states, presumably related to the 



^Now particle data table denote notations /o and ao to 1=0,1 scalars: however particular resonances 

listed there under these names hardly have anything to do with correlators under consideration. 

^ For its charged component (e.g. ud) it is trivial, because two quarks have different flavor: a one-line 

algebra shows why it is so for neutral component as well. 
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famous sigma peak in two pion scattering. As a result, a correlator look quite different 
from what was found in these works. 

Now come better news. The non-locality and difference between a and 5 operators 
are, after all, just lattice artifacts. In continuum limit, both one-loop diagrams become 
identical 

K'[{x,y)^Ki{x,y) (9) 

At T=0 we know that the 5 channel does not have any prominent resonances and is 
dominated by multi-pion states with energy above 1 GeV. And indeed, such features are 



seen in 5 correlator, calculated with Wilson fermions (e.g. [|r^). Thus, we are able to test 
our conjuncture, find qualitative agreement. On this ground, let me speculate that it is 



also true at non-zero T, so results of pO[ do resemble 5 rather than a correlator. If so, 
instead of the well anticipated tt — a degeneracy at they gave actually hint a toward 
more interesting statement, the tt — 5 degeneracy, which implies restoration of both chiral 
symmetries at ('scenario 2'). 

Of course, further studies are needed to clarify this conjuncture. In doing this, it 
would be especially useful to repeat what was done by Gottlieb et al , namely measure 
separately the difference between scalar and pseudoscalar correlators, which has the most 
dramatic temperature (and quark mass) dependence. 

Finally, another way of answering the question we discuss is to study T- dependence 
of the average values of certain four-fermion operators. Selecting those which are related 
to chiral symmetry breaking, one may get rid of perturbative contribution. However, the 
small mass extrapolation is still necessary, as for the quark condensate. 

A quantity sensitive to U{l)a is well known, it is e.g. 't Hooft-type operator 

QtHooft^j.^ =« [iuRUL)idRdL)] » (10) 

(where L,R stand for left and right components of the quark fields). Around or above Tc 
there is no quark condensate, and 'unpaired' instantons have very small density 0{m^f) 
However, in the measurement of quantities like this one, the denominators of the 
quark propagators will take powers of quark masses back, recovering the famous 't Hooft 
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effective interaction. In a sense, tlie operator considered can induce a tunneling event 
by itself, which was not present in the vacuum without it. Unfortunately, it implies in 
practice that small fraction of configurations will produce a large signal: certainly not an 
easy way for measurements. However, this is what one should do, trying to understand 
the question under consideration. 

Let us for completeness also mention some 'alternative SU (Nf) order parameters', 
the expectation values of the following 4-fermion operators 

Oi(r) =« [ul^oUl - dL7odL][L ^ R] » (11) 
02{T) =« [UL^^UL - dLlidL] \L R] » (12) 
OsiT) =« [uLlof'uL - dLi^edL][L ^ R] » (13) 

04(T) =« [ul^^ul - dLirdL][L R] » (14) 

As it is discussed in [^, those enter Weinberg-type sum rules at non-zero temperatures, 
related with the difference between vector and axial correlators. 

5. Can one observe the difference between the two scenarios in 'real' experi- 
ments? 

It is expected that in high energy heavy-ion collisions at RHIC and LHC, the 
system will spend a significant amount of time, of the order of 30 fm/c, in the so called 
'mixed phase' at T Tc. It is not easy to produce quark-gluon plasma, but, after it is 
produced, it also difficult to transform it into ordinary low-density matter. 

A very fascinating possibility, suggested by Bjorken ||23[ and recently studied in 



[0] using 'quenched' scenario, is possible formation of large-amplitude classical pion field 
due to large critical fluctuations. In order to described it one should know which type 
of fields should be included in the effective Lagrangian, describing dynamics of the phase 
transition. Scenario 1 implies that it is a (3-d) variant of sigma model, while scenario 2 
demand the inclusion of (at least) twice more soft modes. 
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I could not figure out any observable signals for both scalars, a and 5 even if 
fluctuations in their directions at the transition is equally strong to that of the pion. 
However, it is clear that scenario 2 implies large fluctuation in rjnon- strange direction, and 
those can in principle be observed. 

First of all, as proven by the WA80 experiment, both vr^,-?] can be extracted in 
the 77 mode, even in large-multiplicity heavy-ion collisions with its huge combinatorial 
background 0. 

Moreover, the rj spectra should have clear signs of the U{1)a restoration, provided 
scenario 2 is the case. In this case t] mesons are easily produced inside the fireball, being 
very light, but then, while trying to leave the system, they should then experience the 
influence of very strong collective potential. Its crude estimate is V ~ (m^ — m^) ~ 
AOOMeV. In contrast to pions, for which it is much smaller, it is comparable to their 
kinetic energy at breakup -Efcmetic ~ 3T ~ AOOMeV. As a result, a significant fraction of 
etas should be trapped inside the fireball, and eventually shifted to much lower kinetic 
energies. This phenomenon, analogous to electrons evaporated from the hot cathode, 
should lead to a peak in the eta spectrum at small pf. 

6. Summary 

Summarizing our discussion, we can repeat that the main question discussed in 
these comments remain open. 

At one hand, it looks plausible that 'phase transition in mathematical sense', the 
singularity, is related with 4 'soft' modes, pions and sigma ones. At another, it becomes 
increasingly clear that physics at T ^ Tc cannot be described in terms of these fields 
alone, and excitations of matter in this region are very different from those at T = 0. 

In particular, there are evidences that parity partners of pions and sigma, 6, rjnon-str 
which are very heavy or even non-existing at T=0, become comparably light. Such ap- 
proximate U{l)a restoration is not however based on 'instanton suppression', but rather 
on new mechanism, a reorganization of the 'liquid' into into the instanton-anti-instanton 

^ Unfortunately, it can only be done statistically, so 'fluctuations in eta direction' on event-per-event 
basis can only be done indirectly, via fluctuations in total number of photons. 
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molecules. 

Finally, if these phenomena take place, large fluctuation of rj production in heavy 
ion collision, especially at small transverse momenta, are predicted. 
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